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Q . Abstract 

^^ ■ We analyze the relation between the Dirac spectrum and the gauge field in SU(3) 

lattice QCD. We focus on how a certain component of the gauge field is related to 
the Dirac spectrum. First, we consider momentum components of the gauge field. 
It turns out that the broad momentum region is relevant for the low-lying Dirac 
spectrum and topological charges. The connection with chiral random matrix theory 
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D \ is also discussed. Second, we consider an SU(2) subgroup component of the SU(3) 
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gauge field. The SU(2) subgroup component behaves like the SU(2) gauge field in 
the low-lying Dirac spectrum. 
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1 Introduction 



5r , The eigenvalue spectrum of tlie Dirac operator is a fundamental quantity in 

QCD. It reflects important properties induced by the gauge field. In particular, 
zero modes of the Dirac operator have special roles on nonperturbative prop- 
erties of QCD. The number of zero modes determines the topological charge 
Q of the background gauge field, as 

nR-nL = NfQ, (1) 

where npt and n^ are the numbers of right-handed and left-handed zero modes, 
respectively [Ij. Also, the existence of zero modes is related to spontaneous 
chiral symmetry breaking, through the Banks-Casher relation. 
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7rp(0) = S, (2) 



where p(0) is the macroscopic spectral density of zero modes and S is the 
absolute value of the chiral condensate [2] . 



Moreover, interestingly, the low-lying Dirac spectrum shows the universal be- 
havior of disordered systems. Such an universal behavior is successfully de- 
scribed by chiral random matrix theory [3lH|5] . The distribution of the low- 
lying Dirac spectrum depends only on the global symmetries of the system, 
and not on the details of the microscopic structure. Once the validity of chi- 
ral random matrix theory is guaranteed, chiral random matrix theory enables 
us to predict the low-lying eigenvalue distribution by a simple and universal 
function. 



It is definitely important to clarify how the gauge field induces these properties 
of the Dirac spectrum. However, it is rather difficult. These properties are be- 
yond the reach of perturbation theory. Although the Dirac equation is written 
in a simple form, the low-energy Dirac spectrum is not easily obtained be- 
cause of the complicated dynamics of the gauge field. In chiral random matrix 
theory, the information about the dynamics of the gauge field is lost. 



The aim of this study is to investigate the relation between the Dirac spectrum 
and the gauge field in lattice QCD. We focus on two kinds of components 
of the gauge field: momentum components and subgroup components. We 
decompose the gauge field into some components, and then calculate the Dirac 
spectrum from the decomposed components. From the change of the Dirac 
spectrum, we analyze how the components of the gauge field are related to 
the Dirac spectrum. Although such a decomposition is gauge dependent, we 
can intuitively understand the role of the gauge field from this analysis. 



In this paper, we calculate the staggered Dirac spectrum in lattice QCD. 
The basics of the staggered Dirac spectrum and the details of the numerical 
simulations are shown in section 2. For the analysis of momentum components, 
we introduce a lattice framework of momentum cutoff in section 3. We discuss 
the resultant Dirac spectrum and the connection with chiral random matrix 
theory. For the analysis of subgroup components, we project the SU(3) gauge 
field onto an SU(2) subgroup component, as demonstrated in section 4. Finally, 
section 5 is devoted to the summary. 



2 The Dirac spectrum in lattice QCD 

2.1 Staggered Dirac operator 

The Euchdean Dirac operator is defined as 

D = Y.hA - ^^7/.^/. (a;)]- (3) 

/^ 

The j-th eigenmode of the Dirac operator is given by 

D4j, = iXj^,. (4) 

The chirality is defined as 

In continuum QCD, the anti-commutation relation, 

{^,75} = 0, (6) 

is satisfied. Because the Dirac operator is anti-Hermitian, A^ is zero or nonzero 
real. Nonzero eigenvalues are always paired as ±«Aj and their chiralities are 
zero. On the other hand, zero eigenvalues are unpaired, and their chiralities 

are +1 or —1. 

The staggered Dirac operator is 



where 77^ (x) is the standard staggered phase factor. For the standard (unim- 
proved) staggered Dirac operator, U^{x) is the SU(3) link variable U^{x). For 
the improved staggered Dirac operator, Ufj_{x) is a combination of the link 
variables. 

The staggered fermion describes the four-fiavor Dirac fermion in continuum 
limit. At finite lattice spacing, however, there is an Ola"^) term which vio- 
lates the fiavor symmetry and Eq. ([6]). Due to the presence of this term, the 
staggered Dirac spectrum does not have exact zero eigenvalues, unlike the con- 
tinuum Dirac spectrum. As a result, the unimproved staggered Dirac operator 



fails to reproduce nontrivial topological sectors PI7|8] . This problem is ap- 
proximately settled by improving the gauge and fermion actions Pll0fllfl2] . 



2.2 Simulation details 



We computed the low-lying eigenvalues of the staggered Dirac operator by 
SU(3) quenched lattice QCD simulations PI7|8]|9]IT0|11|12] . We used the im- 
proved gauge action and the improved staggered Dirac operator. 

For the gauge action, we adopted the 1-loop improved Symanzik gauge ac- 
tion with tadpole improvement. The 1-loop improved Symanzik gauge action 
includes the plaquette term, the 1x2 rectangle term, and the 1x1x1 "par- 
allelogram" term [13]. The tadpole coefficient is defined as the fourth-root of 
the average plaquette value [H] . As listed in Table [H we generated four kinds 
of quenched gauge configurations. The physical volumes of these gauge config- 
urations are approximately the same. The 12^, 16^, and 20^ lattices are used 
for the analysis of eigenvalues and topological charges. The 8^ lattice is used 
for the discussion about chiral random matrix theory. 



For the staggered Dirac operator, we adopted the FATTxASQ operator 
The FAT7xASQ operator is constructed from the FAT7 link variable and the 
ASQ operator [16]. This improved operator strongly suppresses the flavor- 
symmetry breaking lattice artifact. The spatial boundary conditions are peri- 
odic, and the temporal boundary condition is antiperiodic. Because the oper- 
ator is ant i- Hermit ian, the eigenvalues always appear in pairs, ±zAj. We only 
have to calculate the positive eigenvalues \j > 0. It should be understood 
that the negative eigenvalues with the same magnitude exist. We calculated 
the low-lying 30 positive eigenvalues. 

In Fig. [H we show the low-lying 10 positive eigenvalues of one typical gauge 
configuration of /3 = 8.6. The four-fold degeneracy, i.e., the four-fiavor symme- 
try, is satisfied with high precision. In addition, zero eigenvalues are approxi- 
mately reproduced. Such zero modes are called as "would-be" or "near" zero 
modes. Because the negative eigenvalues with the same magnitudes exist, the 

Table 1 

The parameters of quenched gauge configurations. The lattice coupling /3, the 
corresponding lattice spacing a, the lattice volume V, and the configuration number 
A^conf are listed. 



/3 


a (fm) 


V{a^) 


V (fm^) 


A^conf 


9.0 


0.07 


20"^ 


1.44 


50 


8.6 


0.08 


16^ 


1.3^ 


50 


8.3 


0.11 


12^ 


1.3^ 


50 


7.9 


0.16 


8^ 


1.34 


5000 




Fig. 1. The eigenvalue spectrum of the FAT7x ASQ Dirac operator. Only low-lying 

10 positive eigenvalues of one typical gauge configuration are shown. The calculation 

is done on the 16^ lattice with j3 = 8.6. 
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Fig. 2. The scatter plot of the eigenvalue \j and the absolute value of the chirality Xj- 
The low-lying 30 eigenvalues of 50 gauge configurations are plotted. The calculation 
is done on the 16^ lattice with /3 = 8.6. 



topological charge of this gauge configuration is Q = 1. In Fig. O we show the 
scatter plot of all the low-lying 30 eigenvalues of 50 gauge configurations of 
/3 = 8.6. There exists a clear separation between the nonzero modes, which is 
\Xj\ — 0, and the would-be zero modes, which is \xj\ — 0.8. The reason why 
the chirality deviates from ±1 is that it requires a renormalization |17] . 



3 Analysis of momentum components 



In this section, we analyze how the Dirac spectrum is changed when the gauge 
field is restricted to some momentum components. Since the behavior of the 
gauge field depends on the energy scale, different momentum components 
would have different roles. In a different point of view, this analysis relates to 
the connection between the energy scales of the fermion field and the gauge 
field. 
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Fig. 3. The schematic figure of momentum space. 



3.1 Momentum cutoff 



Here, we briefly introduce the lattice framework to remove some momentum 
components of the gauge field [T81IT9] . The framework is formulated as the 
following five steps. 

Step 1. The link variable f/^(x) is generated by Monte Carlo simulation. The 
link variable is fixed with a certain gauge. 

Step 2. The momentum-space link variable U^{p) is obtained by the Fourier 
transformation, as 



U^{p) = —J2Uf,{x) exp{iJ2PuXu) 



where V is the lattice volume. 

Step 3. Some components of U^{p) are removed by introducing a momentum 
cutoff. In the cutoff region, the momentum component is replaced by the free 
gauge field A^ = 0, i.e., the free link variable 



U'rip) = ^ E 1 expitE.P.x,) = 5,0. 



(9) 



For example, as one of the most natural choices, the ultraviolet cutoff Auv is 



defined by 



I (v^>Auv), 



and the infrared cutoff Air is defined by 



v»4''° '^"'"' (11) 



The schematic figure is shown in Fig. |3l 

Step 4. The coordinate-space hnk variable is obtained by the inverse Fourier 
transformation as 



^;w = E^M(p)e^p(-^Ep^^^)- (12) 



Since U'^(x) is not an SU(3) matrix in general, U'^{x) must be projected onto an 
SU(3) element UHx). The projection is realized by maximizing the quantity 



ReTr[{[/^^(x)}tf/;(x)]. (13) 



Step 5. The expectation value of an operator O is computed by using this link 
variable U^{x) instead of Ufi{x), i.e., {0[U^]) instead of {0[U]). 

Using this framework, we analyze how the physical quantity is changed by the 
momentum cutoff. From the resultant change, we can nonperturbatively inves- 
tigate the relation between the physical quantity and momentum components 
of the gauge field. 

While the Dirac spectrum is gauge invariant, a momentum component of the 
gauge field is gauge variant. This is because a gauge transformation is nonlocal 
in momentum space, and some momentum component could become another 
momentum component by a gauge transformation. Thus, we must fix the link 
variable with some gauge in Step 1. Although we can choose any gauge, we 
mainly used the Landau gauge in this paper. 
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Fig. 4. The low-lying Dirac spectrum with the nioinentum cutoff. The original gauge 
configuration is the same as that of Fig. [TJ 
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Fig. 5. The scatter plot of the eigenvalue \j and the absolute value of the chiral- 
ity Xj with the ultraviolet cutoff Auv- The low- lying 30 eigenvalues of 50 gauge 
configurations are plotted. The calculation is done on the 16 lattice with /3 = 8.6. 
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Fig. 6. The scatter plot of the eigenvalue \j and the absolute value of the chirality Xj 
with the infrared cutoff Air. The low- lying 30 eigenvalues of 50 gauge configurations 
are plotted. The calculation is done on the 16^ lattice with /3 = 8.6. 



3.2 Eigenvalue and chirality 



We applied the above framework to the Dirac spectrum. In Fig. HI we show 
the Dirac spectrum with Auv = 3 GeV and Air = 3 GeV. The original gauge 
configuration is the same as that of Fig. [H For Auv = 3 GeV, while all the 
eigenvalues change, the four-fold degeneracy is still satisfied. The would-be 
zero modes which exist in Fig. [1] disappear. This means that the topological 
charge is broken. For Air = 3 GeV, the low-lying 10 eigenvalues become close 
to the eigenvalue of free fermions, which is about 0.39 for this case. Also in 
this case, the topological charge becomes zero. 

In Fig. [5], we show the scatter plot of the eigenvalue and the chirality with 
the ultraviolet cutoff Ayy. As the high- momentum component is removed, the 
eigenvalue spectrum is gradually changed. For Ayy = 6 GeV, the distribution 
is not so different from the original one. For Auv = 3 GeV, all eigenmodes are 
\Xj\ — 0. Therefore, all gauge configurations fall into the trivial topological 
sector. 

The scatter plot with the infrared cutoff Air is shown in Fig. [Hi When the low- 
momentum component is removed, many of the low- lying eigenvalues approach 
the free eigenvalue, which is Xj ~ 0.39 and \xj\ — 0. Compared to the case 
of the ultraviolet cutoff, the changes of the eigenvalues are sensitive to the 
infrared cutoff. This result is consistent with our expectation that the low-lying 
Dirac eigenmode interacts with the low-momentum gauge field. For Air = 3 
GeV, we cannot see would-be zero modes since all eigenvalues are almost 
degenerate. 

To examine discretization effects, we also calculated the 12^ lattice with /3 = 
8.3 and the 20"^ lattice with (3 = 9.0. The results are summarized in Fig. [71 
The distribution of the eigenvalue and the chirality depends on lattice spacing; 
the would-be zero modes are more clearly identified in the case of finer lattice 
spacing. This behavior has been known in the standard lattice QCD pTp2] . 
In contrast, the existence of zero modes and its dependence on the momentum 
cutoff are insensitive to lattice spacing; the would-be zero modes are gradually 
destroyed at Auv = 4 GeV, and then all eigenmodes become nonzero modes 
at Auv = 3 GeV. The lattice regularization is one kind of the ultraviolet 
cutoff, i.e., momentum is limited to less than a half width of the Brillouin zone 
7r/a, as shown in Fig. [31 When lattice spacing is fine enough, the extra-high- 
momentum component near 7r/a gives rise to renormalization effects, but does 
not contribute to physical observables. Actually, this extra-high-momentum 
component does not affect the existence of zero modes. 

When the momentum cutoff is introduced, the clear separation between would- 
be zero modes and nonzero modes is lost. This is because the staggered Dirac 
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Fig. 7. The scatter plot of the eigenvalue \j and the absolute value of the chirality Xj 
"with the ultraviolet cutoff Auv- The results on the 20'* and 12^ lattices are plotted. 

operator does not have exact zero modes. This situation will be improved by 
the overlap operator, which has exact zero modes [201I2T] . 



3.3 Topological charge 



From the number of would-be zero modes, we extracted the topological charge 
Q of the background gauge field. The topological charge is given through 
Eq. ([1]). We set the criterion for would-be zero modes as \xj\ > 0.6, and for 
nonzero modes as \xj\ < 0.6. This criterion clearly separates would-be zero 
modes in the original case, as shown in Fig. [2J For the cases in which the 
clear separation between would-be zero modes and nonzero modes is lost, this 
criterion could be subtle. The results are summarized in Table |2l 

For a consistency check, we also extracted the topological charge in another 
way. We directly calculated the topological charge by 






(14) 



Its discretized form is 



Q 



±4 

E E 






Tr[f/^^(x)f/^p(x)], 



(15) 



10 



Table 2 

Topological charge Q estimated by would-be 
by Eq. (|15p is shown in the parenthesis. The 
are shown, 



zero modes. The geometrical estimate 
results of the 50 gauge configurations 



V{a^) 


Data 


IQI=o 


IQI = i 


IQI = 2 


IQI=3 


16^ 


No Cut 


28 (29) 


18 (17) 


4(4) 


0(0) 




Auv = 6 GeV 


35 (27) 


12 (17) 


3(6) 


0(0) 




Auv = 4 GeV 


50 (44) 


0(6) 


0(0) 


0(0) 




Auv = 3 GeV 


50 (50) 


0(0) 


0(0) 


0(0) 




Air = 1 GeV 


29 (27) 


17 (19) 


4(4) 


0(0) 




Air = 2 GeV 


32 (43) 


13(7) 


5(0) 


0(0) 




Air = 3 GeV 


49 (50) 


1(0) 


0(0) 


0(0) 


20^ 


No Cut 


28 (28) 


22 (22) 


0(0) 


0(0) 




Auv = 4 GeV 


50 (50) 


0(0) 


0(0) 


0(0) 




Auv = 3 GeV 


50 (50) 


0(0) 


0(0) 


0(0) 


12^ 


No Cut 


28 (27) 


20 (19) 


2(2) 


0(0) 




Auv = 4 GeV 


49 (44) 


1(6) 


0(0) 


0(0) 




Auv = 3 GeV 


50 (50) 


0(0) 


0(0) 


0(0) 



Table 3 

Topological charge Q in the Coulomb gauge. The results of the 50 gauge configu- 
rations of the 16^ lattice are shown. 



V (a^ 



Data 



|Q|=0 \Q\ = l \Q\=2 |Q|=3 



16^ 



No Cut 


28 (29) 


18 (17) 


4(4) 


0(0) 


Air = 1 GeV 


29 (27) 


17 (20) 


4(3) 


0(0) 


Air = 2 GeV 


35 (31) 


13 (14) 


1(5) 


1(0) 


Air = 3 GeV 


49 (50) 


1(0) 


0(0) 


0(0) 



where U^u{x) is the plaquette. We measured this topological charge after the 
APE smearing steps [8|22] . The smearing coefficient and the number of the 
smearing step are set to 7.0 and 200, respectively. The results are shown in the 
parentheses in Table [2l The obtained topological charge is roughly consistent 
with the estimate by the would-be zero modes, especially in the original lattice 
QCD ("No Cut"), as expected. These two estimates deviate from each other 
in some cases with the momentum cutoff, e.g., Auv = 4 GeV. This is partly 
because it is difficult to identify the would-be zero modes in these cases, and 
partly because the APE smearing steps change the local topological structure. 

We can see that the number of the topological charges is gradually reduced in 
the both cases of the ultraviolet and infrared cutoffs. The topological charge 
is almost completely destroyed at Air = 3 GeV or Auv = 3 GeV. This means 
that the topological charge can be destroyed by removing the momentum 
region to some extent, regardless of whether it is ultraviolet or infrared. In 
other words, there is no typical energy scale of the gauge field inducing the 
topological charge. This is different from the case of color confinement, which 
is induced only by the narrow low-momentum component of the gauge field 



11 






Fig. 8. The probability distribution of the lowest eigenvalue Amin in the Q = 
sector. The histograms are the lattice data of the 8^ lattice with /3 = 7.9. The solid 
curves for No Cut and Auv = 4 GeV are the best-fit functions of Eq. (J16p from 
chiral random matrix theory. 

We can study any other gauges in the same manner, instead of the Landau 
gauge. In Table [3l we show the numerical results with the Coulomb gauge. 
The results are almost the same as those of the Landau gauge. 



3.4 Connection with chiral random matrix theory 



Next, we compare the lattice data with a prediction of chiral random matrix 
theory, using the gauge configurations oi (3 = 7.9 and V = 8"^. This lattice 
volume is large enough to discuss the connection with chiral random matrix 
theory, i.e., in the so-called e- regime [11]. Here, we analyze the Q = sector 
with the criterion for nonzero modes as \xj\ < 0.1. 

We calculated the probability distribution of the lowest nonzero eigenvalue 
Amin- In chiral random matrix theory, the probability distribution of the lowest 
eigenvalue is given as 



P{K 



-22. /4 

mm' 



'^min — ^miri'^ ^^ 



(16) 
(17) 



for the SU(3) quenched gauge field with Q = [21]. If one directly calculates 
the chiral condensate in infinite volume, this function is parameter free. In 
this study, we treated S as a fitting parameter. The lattice data and the best- 
fit result of Eq. ( IT6l) are shown in Fig. |H1 In the original lattice QCD (No 
Cut), the lattice data is well reproduced by chiral random matrix theory. The 
best-fit value is S ~ 0.0022 x a'^. 



For Auv = 4 GeV, the qualitative behavior does not change. The lattice 
data is still reproduced by Eq. (TT6|l . Thus, the high- momentum component 
of the gauge field is irrelevant for the universality of the lowest eigenvalue 
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distribution. The best-fit value of E is clianged by tlie ultraviolet cutoff. This is 
because the chiral condensate is renormalization-group variant and ultraviolet 
divergent. By definition, the ultraviolet divergent quantity depends on the 
ultraviolet cutoff, although this might be irrelevant for phenomenology. 

For Air = 2 GeV, the lowest eigenvalue distribution becomes to concentrate 
on the vicinity of the free eigenvalue. The lattice data drastically deviates 
from Eq. f ll6p . Therefore, the low- momentum component of the gauge field is 
crucial for the validity of chiral random matrix theory. In other words, the low- 
momentum component induces the strong-interacting and disordered nature 
of the gauge field. 

In a previous work, the chiral condensate was directly analyzed by the lattice 
framework of the momentum cutoff p5] . It suggests that the zero-momentum 
component of the gauge field has a large contribution to the chiral condensate. 
To check it in the present study, we calculated the Dirac spectrum with the 
small infrared cutoff Air ~ 0.1 GeV. For Air ~ 0.1 GeV, the lattice data can 
be fitted by Eq. fITB]) with the best-fit value E ~ 0.0019 x a~^. The contribution 
of the zero-momentum component is indeed large, about 15% of the total. This 
is qualitatively consistent with the previous work. Note, however, that this is 
quantitatively different from the previous work because of many systematic 
differences. 



4 Analysis of subgroup components 



Although there are many similarities between the SU(2) and SU(3) gauge 
theories, there are several differences because of an additional symmetry of 
SU(2). An SU(3) group element includes SU(2) subgroup elements. In this 
section, we consider whether an SU(2) subgroup component of the SU(3) gauge 
field shows the same behavior as the SU(2) gauge field. 



4.1 Subgroup projection 



The framework of the subgroup projection is formulated as below. 

Step 1. The link variable U^{x) is generated by Monte Carlo simulation. The 
link variable is fixed with a certain gauge. 

Step 2. From the SU(3) link variable. 
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Table 4 

Topological charge Q estimated by would-be zero modes. The geometrical estimate 
by Eq. (|15p is shown in the parenthesis. The results of the 50 gauge configurations 
with /3 = 8.6 are shown. 



V ia'' 



Data 



|Q|=0 \Q\ = l \Q\ = 2 |Q|=3 



16^ SU(3) 28 (29) 18 (17) 4 (4) (0) 

SU(2) projected 38 (35) 10 (11) 2 (3) (1) 



'^t/n f/i2f/i3^ 



U,{x) 



U21 U22 U23 
U31 f/32 U33 



the SU(2)-projected link variable is generated from its subgroup element, as 



Uf,{x) 



^UnUi20^ 



U21 U22 
1 



(19) 



Step 3. The expectation value of an operator O is computed by using this link 
variable u^{x) instead of U^{x), i.e., {0[u]) instead of {0[U]). 



4.2 Numerical results 



The numerical analysis was performed in the same way as the previous section. 
Also for the subgroup projection, we used the Landau gauge. We can expect 
that other possible embeddings of the SU(2) subgroup give rise to the same 
results because the matrix elements are uniform in the Landau gauge. 

We show the scatter plot of the eigenvalues and the chiralities in Fig. [9l and 
the list of the topological charges in Table HI As for would-be zero modes and 
topological charges, the SU(2)-projected gauge field is qualitatively similar to 
the original SU(3) gauge field. The Dirac spectrum has would-be zero modes, 
and the SU(2) topological charge exists. Note that the number of would-be 
zero modes in each gauge configuration is different from that of the original 
SU(3) theory. 

The low-lying Dirac spectrum depends on the symmetry of the fermion field 
[26] . The eigenvalue distributions are completely different between SU(3) and 
SU(2). For the staggered Dirac operator with the SU(2) quenched gauge field, 
chiral random matrix theory predicts that the lowest eigenvalue distribution 
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Fig. 9. The scatter plot of the eigenvalue Xj and the absolute value of the chirality 
Xj of the SU(2)-projected gauge field. The low-lying 20 eigenvalues of 50 gauge 
configurations are plotted. The calculation is done on the 16^ lattice with /3 = 8.6. 
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Fig. 10. The lowest eigenvalue distribution of the SU(2)-projected gauge field in the 
Q = sector. The solid and dashed curves are Eqs. (pOj) and p^ . respectively. 



IS 



P{K 



■^min COSll ^min 



sinh Zr, 



./2 



(20) 



in the Q = sector [24]. This distribution is indeed obtained in SU(2) lattice 
QCD [27]. In Fig. [TOl we show the lattice data of the SU(2)-projected gauge 
field. The solid and dashed curves are the predictions of chiral random matrix 
theory for SU(2) and SU(3), i.e., Eqs. f l20|) and f lT6|l . respectively. The lattice 
data is reproduced by the prediction for SU(2), not for SU(3). Therefore, 
the SU(2) subgroup component of the SU(3) gauge theory shows the same 
universal behavior as the SU(2) gauge theory. The best-fit value of the chiral 
condensate is E ~ 0.0013 x a~^. The ratio of this SU(2)-projected chiral 
condensate to the original SU(3) chiral condensate is 0.6, which is roughly 
understood as the ratio of the numbers of color. 
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5 Summary 



We investigated the relation between the Dirac spectrum and components of 
the gauge field. We considered momentum components and a subgroup com- 
ponent of the gauge group. As for momentum components of the gauge field, 
the Dirac zero mode and the topological charge are gradually destroyed by 
the ultraviolet or infrared cutoff. The low-momentum component is especially 
crucial for the validity of chiral random matrix theory. As for an SU(2) sub- 
group component of the SU(3) gauge group, it behaves like the SU(2) gauge 
field. 

The relation between zero modes and the gauge field is relevant for the topo- 
logical structure of QCD vacuum, e.g., the topological charge density distribu- 
tion. It would be instructive to perform the similar analysis with the overlap 
operator. 
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